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Abstract 

We propose a new model of modified F{R) gravity theory with 
the function F{R) = (1//3) arcsin(/3i?). Constant curvature solutions 
corresponding to the flat and de Sitter spacetime are obtained. The 
Jordan and Einstein frames are considered; the potential and the mass 
of the scalar degree of freedom are found. We show that the flat space- 
time is stable and the de Sitter spacetime is unstable. The slow-roll 
parameters e, g, and the e-fold number of the model are evaluated in 
the Einstein frame. The index of the scalar spectrum power-law Ug 
and the tensor-to-scalar ratio r are calculated. Critical points of au¬ 
tonomous equations for the de Sitter phase and the matter dominated 
epoch are found and studied. We obtain the approximate solution of 
equations of motion which is the deviation from the de Sitter phase 
in the Jordan frame. It is demonstrated that the model passes the 
matter stability test. 


1 Introduction 

It is possible to describe the inflation and the present time universe accel¬ 
eration if one modihes the Einstein-Hilbert (EH) action of general relativity 
(GR). We propose the particular model of the F{R) gravity with the help 
of replacing the Ricci scalar by the function F{R) = (1//9) arcsin(/5R) in EH 
action, where [3 is the parameter with the dimension of (length)^. Thus, we 
introduce the fundamental length \/^ which goes probably from quantum 
gravity. The F{R) gravity models may describe the evolution of the universe 
without introducing Dark Energy (DE) (Appleby, Battye and Starobinsky 
2010, Capozziello and Faraoni 2011, Nojiri and Odintsov 2011). In such 
models the cosmic acceleration occurs due to modihed gravity. Therefore, 
F{R) gravity models can be an alternative to A-Cold Dark Matter (ACDM) 
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model as new gravitational physics is considered. The ACDM model has 
a problem with the explanation of the smallness of the cosmological con¬ 
stant A. It should be mentioned that the form of the function F{R) has to 
be derived from the fundamental theory (string, M-theory) which is absent. 
Therefore, different F{R) gravity models that satishes the general conditions 
are of interest. The motivation for this work is to consider new F(R) model 
which meets requirements such as quantum and classical stabilities, it passes 
matter stability test and describes inflation of universe etc. 

It should be noted that the hrst successful models of F{R) gravity were 
given in Starobinsky 1980, Hu and Sawicki 2007, Appleby and Battye 2007, 
Starobinsky 2007, Nojiri and Odintsov 2007a, Nojiri and Odintsov 2008, 
Cognola et al 2008. Some F{R) gravity models were introduced in Deser 
and Gibbons 1998, Capozziello and Faraoni 2011, Kruglov 2013, Kruglov 
2014a, Kruglov 2014b and in other publications. F{R) gravity models are 
phenomenological models that may describe different eras and the evolution 
of the universe. The hrst F{R) gravity model was introduced in Starobinsky 
1980 that gives the self-consistent description of the inhation. 

The paper is organized as follows. We formulate the model with one 
dimensional parameter in section 2. It is shown that the classical and 
quantum stabilities take place in the model under consideration. We obtain 
the constant curvature solutions corresponding to hat spacetime, Rq = 0, 
and to the Schwarzschild-de Sitter spacetime, /3Ro ~ 0.919. In section 3 the 
scalar-tensor formulation of the model is investigated (in the Einstein frame). 
The potential and the mass of the scalar degree of freedom (scalaron) are 
found. We obtain the slow-roll parameters e, rj, and the e-fold number of 
the model in section 4. The index of the scalar spectrum power-law Ug and 
the tensor-to-scalar ratio r are calculated. Critical points of autonomous 
equations for the de Sitter phase and the matter dominated epoch are found 
in section 5. The approximate solution of equations of motion in the Jordan 
frame corresponding to the deviation from the de Sitter phase is obtained in 
section 6. We show in section 7 that the model passes the matter stability 
test. Section 8 is devoted to a conclusion. 

The Minkowski metric ? 7 ^,^=diag(-l, 1, 1, 1) is used and c=h=l is as¬ 
sumed. 


2 



2 The Model 

Let us consider a new model of arcsin-gravity with the Lagrangian density 


C 


2^2 


F{R) 


1 


1 


arcsin(/9i?) , 


( 1 ) 


where k = Mpl, Mpi is the reduced Planck mass, /9 has the dimension 
of (length)2, and the action without matter is given hy S = f (Rx^/^C. 
At /Si? 1, we have arcsin(/Si?) /Si?, and we arrive at the EH action. 
The equation E(0) = 0 holds, corresponding to the flat space-time without 
cosmological constant. GR passes local tests and we imply that at the present 
time the low curvature regime occurs, /Si? 1. We will describe the inflation 
and universe evolution in the model suggested. For the classical stability the 
inequality F'{R) > 0 (the prime means the derivative with the respect to 
the argument) is required (Appleby, Battye and Starobinsky 2010) which is 
satished if /Si? < 1, 


F'(R) 


1 

\/i - (mr 


> 0 . 


( 2 ) 


Quantum stability claims the inequality F'\R) > 0 (Appleby, Battye and 
Starobinsky 2010), that becomes in our model as follows: 


F"(i?) = 


/S^i? 


[1 - (/ Si ?) 2]"/2 

and it is also satished at 0 < /Si? < 1. 


> 0 , 


( 3 ) 


2.1 Constant Curvature Solutions 

If the Ricci scalar i? is a constant, i? = i?o, equations of motion (Barrow, 
1983) become 

2F(i?o) = i?oR'(i?o), (4) 

and are given, in the model with the Lagrangian density (1), by 

2-y/l — (/Si?o)2 arcsin(/Si?o) = /Si?o. (5) 

We note that constant curvature solutions correspond to the extremum of 
the effective potential. It should be mentioned that the general conditions 
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for multiply de Sitter solutions in F{R) gravity were discussed in Cognola et 
al 2005, and Cognola et al 2009. The fact of the appearance of such solutions 
is known for F(R) gravity theories with many parameters. But the model (1) 
has only one parameter P that is, in our opinion, an attractive feature of the 
model. Therefore, it is of interest to analyze one more F{R) gravity model. 

Eq. (5) possesses two solution, i?o = 0 corresponding to the flat space- 
time, and non-trivial solution /SRq ~ 0.919. We will show that the last solu¬ 
tion goes with the Schwarzschild-de Sitter spacetime and with the maximum 
of the effective potential in the Einstein frame. The constant curvature solu¬ 
tions describe the acceleration phase which is future stable if the inequality 
F'(Rq)/F"(Rq) > Rq occurs (Muller, Schmidt and Starobinsky 1988), and 
we have from Eqs. (2),(3) 


l-{/3Rf>WR)\ ( 6 ) 

which is equivalent to f3R < l/\/2 ^ 0.707. Thus, the solution i?o = 0 obeys 
Eq. (6) and the flat spacetime is stable. The second constant curvature 
solution /3Ro ~ 0.919 does not satisfy Eq. (6), leads to unstable de Sitter 
spacetime, and describes the inflation. 


3 The Scalar-Tensor Formulation 


Now we investigate the model in Einstein’s frame performing the conformal 
transformation of the metric (Magnano and Sokolowski 1994) 

g^u = F {R)gfj,i, = . ==g^iy. (7) 

^1 - [fjRy 

Then the Lagrangian density in Einstein’s frame becomes 


£ = ^ - V(4,). (8) 

Here the Ricci scalar R in the Einstein frame is calculated in new metric (7) 
and the scalar held 0 is 


4>{R) 


\P1k 


InF'(R) 


\P1k 




(9) 
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Figure 1: The function ncj) versus I3R. 


The plot of the function ncj) is presented in Fig. [T] The potential V{(j)) is 
given by 


V{R) 


RF'{R) - F{R) 
2k‘^F'^{R) 


1 

2/3k^ 


l3R\Jl — {13Ry — (^1 — {(3rY^ arcsin(/9i?) 


( 10 ) 


The plot of the function versus 13R is given in Fig. |2] and the plot 

of function f3iR‘V versus is represented by Fig. |3l The extremum 
of the potential, V'{R) = 0, with the help of Eq. (10) leads to Eq. (4). 
The potential (10) possesses the minimum at i? = 0 and the maximum at 
/3i?o ~ 0.919. The flat space-time {R = 0) is the stable state and the state 
with the curvature Rq ~ 0.919//9 is unstable. 

We obtain the mass squared of a scalaron (scalar degree of freedom) from 
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Figure 2: The function I3k^V versus (3R. 


Eq. (10), 



mR) \ 

F'^{R)) 


X arcsina; 


( 11 ) 


where x = I3R. The plot of the function versus x = I3R is given by 
Fig. m One can verify that < 0 for the constant curvature solution 
Rq ~ 0.919//3, and, therefore, this solution corresponds to unstable state as 
it was mentioned before. It follows from Eq. (11) that at 0.529 < < 1 


m 


we have non-stable states, < 0. The stability of the de Sitter solution 
F(R) gravity models was hrst studied by Muller, Schmidt and Starobinsky 
1988. To pass the Solar system tests the value should be positive and 
big. If the value I3R is small the mass rritj, is big according to Fig. |4] and 
corrections to the Newton law are negligible. 
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Figure 3: The function versus K(f). 


To assure that corrections of F{R) gravity model are small as compared 
to GR for i? 3> -Ri, where Ri is a curvature at the present time, the relations 

\F{R)-R\<R, |F'(R)-1|<1, |RF"(R)|<1 (12) 

should hold (Appleby and Battye 2010). As arcsinx >a;atl>a;>0 
(x = fiR), the hrst inequality in Eq. (12) becomes arcsinx < 2x, and it is 
satished at 1 > x > 0. The second inequality in Eq. (12) is equivalent to 
X < \/3/2 0.866 (as F\R) > 1 for 0 < x < 1). The third inequality in 

Eq. (12) holds at 0 < x < 0.655. As a result, all Eqs. (12) are satished if 
0 < X < 0.655. 
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Figure 4: The function versus (3R. 


4 Slow-Roll Cosmological Parameters 

The slow-roll parameters are given by (Liddle and Lyth 2000) 


(m' 


h(0) = 


V"(^) 

' V(c/>) ■ 


(13) 


When conditions \ri{(j)) \ <C 1, e(0) 1 hold the slow-roll approximation takes 

place. From Eqs. (10),(11) we hnd the slow-roll parameters as follows: 


1 

e = - 
3 


RF\R)-2F{R) 
RF'{R) - F{R) 


1 


arcsm x — x 


\/l — 


x^ arcsm x — x 


( 14 ) 












( 15 ) 


F'\R) + F"{R) [RF\R) - AF{R)] 

F"{R) [RF'{R) - F{R)] 

2 — 4a;\/r^-^arcsina;^ 

3a; — \/l — arcsin x^ 

The plots of the functions e, r] are given in Fig. El Fig. [HI The equation e = 1 




Figure 5: The function e versus I3R. 

has the solution x ^ 0.766. It follows from Fig. |5]that at 1 > fiR > 0.766 
the inequality e < 1 holds. The equation \r]\ = 1 is satished at x ~ 0.516, 
X ^ 0.544 and x ^ 0.925. At 0.544 > (3R > 0.516 and at 1 > f3R > 0.925, 
we have the result \r]\ < 1. Therefore, the slow-roll approximation, e < 1 and 
|? 7 | < 1, takes place at 1 > (]R > 0.925. 

The age of the inflation can be obtained by calculating the e-fold number 
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Figure 6; The function rj versus I3R. 


(Liddle and Lyth 2000) 


1 r<t> VU) 

We find, from Eqs. (9),(10), the number of e-foldings 


r^o X 1 


— x"^ arcsin x — x'^ 

1 dx 

^ ^end — 

a;2) 1 

^2^/l — x"^ arcsin X — xj 


(16) 


(17) 


were Xend = fdRend corresponds to the time of the end of inflation when e or 
|? 7 | are close to 1. Thus, inflation ends when slow-roll conditions are violated. 
We obtain the amount of inflation W ~ 9.7 at Xq = 0.9999 and Xend = 0.92, 
and, therefore, the model can describe the inflation. It should be noted that 
it is required around 60 e-foldings of inflation to solve the flatness and horizon 
problems. 
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Due to density perturbations the index of the scalar spectrum power-law 
is given by the relation (Liddle and Lyth 2000) 

= 1 — 6e -I- 27 ]. (18) 

Using Eqs. (14),(15), the plot of the function of rig versus (3R is represented 
in the Fig.7. The tensor-to-scalar ratio is dehned by Liddle and Lyth 2000, 



Figure 7: The function Ug versus /3R. 

r = 16e. The PLANCK experiment gives the result (Ade et al 2014) 

n, = 0.9603 ± 0.0073, r < 0.11. (19) 

One can see from Fig. 7 that the experimental value of Ug is not satished. 
But the bound for tensor-to-scalar ratio r < 0.11 is satished for 0.933 > 
(3R > 0.905. As a result, the model suggested can give only approximate 
description of cosmology in the Finstein frame. 
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It was stated in Bamba at al 2014, Bamba and Odintsov 2015 that cos¬ 
mology in the Einstein and Jordan frames can be different and these frames 
are physically non-equivalent. As a result, theories in the Einstein and Jor¬ 
dan frames may be considered as different cosmological theories. One can 
recalculate inflationary parameters in the F{R) frame (see Bamba at al 2014). 


5 Critical Points of Autonomous Equations 

To investigate critical points of equations of motion in the Jordan frame, it 
is useful to introduce the dimensionless parameters (Amendola et al 2007) 
which become 

F'{R) XX 


Xi = 


X2 = 


HF'{R) H{1 
F{R) 


x^ 


x^ arcsm x 


6 F'{R)H‘^ 


613m 


m = 


H „ 

^3 - + 2 , 

RF''{R) x'^ 


( 20 ) 


r = 


F'{R) 
RF'{R) _ X 3 
' F{R) ~ 1^2 ' 


(1 


x^ 


X 


vT 


x^ arcsm x 


( 21 ) 


where if is a Hubble parameter, x = (3R, and the dot over the variables means 
the derivative with respect to the time. The deceleration parameter q is given 
by g = 1 — X 3 . Equations of motion in the absence of the radiation, prad = 0, 
with the help of Eqs. (20), (21) can be written in the form of autonomous 
equations (Amendola et al 2007). One can investigate the critical points of 
the system of equations by the study of the function m{r) which shows the 
deviation from the ACDM model. The plot of the function m(r) is presented 
by Fig. [HI The de Sitter point Pi (Amendola et al 2007), in the absence of 
radiation, X 4 = 0 , corresponds to the parameters xi = 0 , X 2 = — 1 , X 3 = 2 
[H = 0, m = i?/12, r = —2). The point Pi corresponds to the constant 
curvature solutions that may be verihed using Eqs. (5),(20). The effective 
equation of state (EoS) parameter, Wg//, and the parameter of matter energy 
fraction, are given for this point by 


Weff = -1 - 2H/ (3ii^) = -1, = 1 - 3)1 - 2)2 - 3)3 = 0, 


( 22 ) 
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Figure 8 : The function m{r) 


which correspond to de Sitter phase. This point mimics a cosmological 
constant and the deceleration parameter becomes q = —1. The constant 
curvature solution x ~ 0.919 corresponds to unstable de Sitter space as 
1 < m{r = — 2 ) 5.4 (Amendola et al 2007). 

For the critical point P 5 ( 0:3 = 1/2), m ~ 0, r —1, and EoS of a matter 
era is tCe// = 0 (a = Then we have a viable matter dominated epoch 

prior to the late-time acceleration (Amendola et al 2007). The equation 
m = —r — 1 has the solution m = 0 , r = — 1 , P = 0 , corresponding to 
the point P 5 . One can verify with the help of Eq. (21) (see Fig. [H]) that 
m'(r = — 1 ) = 0 . As a result, the condition m'(r = — 1 ) > —1 holds and 
we have the standard matter era (Amendola et al 2007). Therefore, the 
correct description of the standard matter era occurs in the model under 
consideration. To investigate the possibility of late-time acceleration (DE) 
in the model, one needs to solve and analyze autonomous equations. The 
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unification of inflation with DE has been proposed hrst in F{R) gravity by 
Nojiri and Odintsov 2003, and Nojiri and Odintsov 2007b. 


6 Equations of Motion in the Jordan Frame 
and Their Approximate Solutions 


It should be mentioned that the change of the frame includes the change 
of the time scale and, therefore, the values of the slow-roll parameters are 
different in the frames. Now we consider equations of motion in the Jordan 
frame and their approximate solutions. If one adds to (1) the Lagrangian of 
the matter with the energy-momentum tensor we obtain equations of 

motion 

R,.F'{R) - ^g,.F{R) + pF\R) - V,V,F\R) = (23) 

where is a covariant derivative. We consider the homogeneous, isotropic 
and spatially flat Friedmann-Robert son-Walker (FRW) cosmology with the 
line element 

ds^ = —dt^ -\- a^{t) {dx^ + dy"^ + dz^'^ . (24) 

Taking the trace of the left and right sides of Eq. (23) one hnds 

RF\R) - 2F{R) + ‘ig'^NaVyF\R) = k^T, (25) 

where T = g^'^Tj;^^. For the FRW metric (24) we obtain the expression 


g^Na^pR = {-g)- 

Then with the help of Eqs. 
curvature 


1/2 


d,. 




= -R- 3HR. 


(26) 


(1), (26) one hnds the equation for the scalar 


R 


^arcsin(/3R)-3f^ + 3Jf4 


/S 


dt 


/i - mr 


K^r. ( 27 ) 


The Ricci scalar can be expressed through the Hubble parameter as follows: 


R = 12H^ + 6H. 


(28) 


We consider the case T = 0, corresponding to the absence of the matter 
or the traceless of the electromagnetic helds. For the case of the constant 
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curvature Eq. (27) converts to Eq. (4). There exists the exact solution to 
Eqs. (27), (28) 

Ro = 0, ^ (29) 

which is the same as the solution in GR. The solution (29) with the flat 
spacetime corresponds to the radiation era. To describe all eras one should 
solve the system of nonlinear equations (27), (28). We consider the approx¬ 
imate solution to this system of equations which is the deviation from the 
constant curvature solution /3Ro ~ 0.919. Thus, we expand R = Rq + Ri that 
is the deviation from de Sitter phase, and assuming i?i i?o- Linearizing 
Eq. (27), we obtain 


2(/5Ro)' - 1 


Ri — 3/3‘^Rq (^Ri -\- SHqRi^ — 0, 


(30) 


where Hq = Rq/12. Solutions to the linear Eq. (30) is in the form Ri = 
Rexp(n±t) with 


-9i7o/jRo ± ^8177g(/jRo)2 + 12R,[2{I3R,Y - 1] 

6f3Ro 


(31) 


The physical solution with the decreasing curvature is Ri = Aexp{nR) which 
for /3Ro ~ 0.919 becomes 


Ri = A exp 


/ 1.066t\ 

[~^r 


(32) 


The linear equation (30) does not fix the amplitude A that is small compared 
to Rq so that /3A 0.919. Eq. (32) shows the rate of decreasing the 

curvature in the de Sitter point. Thus, the approximate solution to Eq. (27) 
at T = 0 for small deviation from the de Sitter phase is I3R = 0.919 -|- 
/SRexp 1.066f/\/(d^. Thus, we propose F{R) gravity model describing 
universe inflation and at weak curvature it becomes GR. The model under 
consideration is eternal F(R) inflation model. Some models of such sort were 
discussed by Nojiri and Odintsov 2015. To describe all phases of the universe 
evolution one needs to solve the system of equations (27), (28) analytically 
or numerically. We leave this for further investigations. 
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7 Matter Stability 


We follow the method of Dolgov and Kawasaki 2003 to investigate the matter 
stability. For weak gravity objects the Minkowski metric (flat) can be used 
and the approximate relation 5 '“'^VqV ^3 ~ holds. When R is uniform 

(for spatially constant distribution) equation (25) becomes 

- 3F(2) {R)R - 3F(3) {R)R^ + {R)R - 2F{R) = k^T, (33) 

where = (FF{R)/dR^. Let us consider a perturbative solution R = 

Rq + Ri with Ri being the perturbed part, |i?i| |.Ro|)- According to GR 

the curvature in the lowest order is given by Rq = —k^F inside the matter 
and i?o = 0 outside the matter. Following to Nojiri and Odintsov 2003, 2011, 
from equation (33), one hnds 


R0 + -R1+ I (Rq + 2 RqR] 


2 F(Ro 


Rn 


1 + FF[r^) 


3F(2)(Ro) 


where 


G(Ro) 


_p(3)2 _ _p(2) p{i) 
_P(2)2 



G(Ro)i?i, 

(34) 


(RqF^"^^ - FF^ F(2) ^ (^2F - RoFF - 
3_P(2)2 


(35) 


The matter is unstable if U{Rq) > 0 as Ri exponentially increases in the 
time. From Eq. (1) we obtain the function U{Ro): 




(1 + 2x'^)\/l — 2 x^ — 1 

^ X ' 

(36) 

where x = /3Ro. One can check that V{Ro) < 0 for 0 < a; < 1. Thus, for 
almost constant curvature, Ro ~ 0, there is no matter instability, U{Ro) < 0. 
As a result, this indicates on stability of the gravitational system and the 
model passes the matter stability test. 


V{Ro) = 


3/3 


2 arcsin(a;) — 


X 


VT 


X 


x^ 
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8 Conclusion 


We propose a new model of modified F{R) gravity representing the effective 
gravity model which can describe the evolution of the universe. The constant 
curvature solutions, (3Ro = 0, (3Ro = 0.919, were obtained that correspond 
to the flat spacetime and the de Sitter spacetime, correspondingly. The de 
Sitter spacetime gives the acceleration of the universe and corresponds to the 
inflation. The flat space-time is stable but the de Sitter spacetime is unstable 
in the model and it goes with the maximum of the effective potential in the 
Einstein frame. The Jordan and Einstein frames were considered and we have 
obtained the potential and the mass of the scalar degree of freedom. The 
slow-roll parameters e, r] and the e-fold number of the model were evaluated. 
The model gives e-fold number Ng ^ 9.7, characterizing the age of inflation, 
in the Einstein frame. We note that the values of the slow-roll parameters 
and e-foldings are different in the Einstein and Jordan frames. We show 
by the analysis of critical points of autonomous equations that the standard 
matter era exists and the standard matter era conditions are satisfied. We 
have obtained the approximate solution of equations of motion which is the 
deviation from the de Sitter phase. To verify that the model can be consistent 
with the accelerating expansion of the present universe if the curvature is 
small, one should solve the system of equation (27), (28). We leave such 
investigation for the future. It should be noted that the effective Newton 
constant in F{R) gravity models is given by Ge// = G/F'{R) and in the 
model (1) becomes Ggff = G^l — {(3RY- Thus, in the small curvature 
regime, /Si? 1, the corrections to the Newton law are negligible. It was 
shown that the model passes the matter stability test. The model may be 
alternative to GR, and can describe early-time inflation. The possible future 
singularities in this model were not investigated that we leave for further 
study. 


References 

[1] Ade, P. A. R., et al 2014, Astron. Astrophys., 571, A16 
(arXiv:1303.5083), (arXiv:1303.5082), (arXiv:1303.5076). 

[2] Amendola, L., et al. 2007, Phys. Rev. D, 75, 083504 

(arXiv:gr-qc/0612180). 


17 



[3] Appleby, S. A., and Battye, R. A. 2007, Phys. Lett. B, 654, 7 
(arXiv:0705.3199). 

[4] Appleby, S. A., Battye, R. A., and Starobinsky, A. A. 2010, JCAP, 
1006, 005 (arXiv:0909.1737). 

[5] Bamba, K., et al, 2014, Phys. Rev. D 90, 124061 (arXiv: 1410.3993). 

[6] Bamba, K., Odintsov, S. D. 2015, Symmetry, 7, 220 (arXiv:1503.00442). 

[7] Barrow, J. D., and Ottewill, A. C. 1983, J. Phys. A, 16, 2757 

[8] Capozziello, S., and Faraoni, V., 2011, ’’Beyond Einstein Gravity: A 
Snrvey of Gravitational Theories for Gosmology and Astrophysics”, 
Springer Science+Bnsiness Media B.V.: New York 

[9] Gognola, G., et al 2005, JGAP, 010, 0502 (arXiv:hep-th/0501096). 

[10] Gognola, G., et al 2008, Phys. Rev. D 77, 046009 (arXiv:0712.4017). 

[11] Gognola, G., et al 2009, Phys. Rev. D 79, 044001 (arXiv:0810.4989). 

[12] Deser, S., and Gibbons, G. W. 1998, Glass. Qnant. Grav., 15, L35 
(arXiv:hep-th/9803049). 

[13] Dolgov, A. D., and Kawasaki, M. 2003, Phys. Lett. B 573, 1 
(arXiv:astro-phys/0307285). 

[14] Hn, W., and Sawicki, 1. 2007, Phys. Rev. D, 76, 064004 
(arXiv:0705.1158). 

[15] Krnglov, S. 1. 2013, Int. J. Theor. Phys., 52, 2477 (arXiv:1202.4807). 

[16] Krnglov, S. 1. 2013, Int. J. Mod. Phys. A, 28, 1350119 (arXiv:1204.6709 

)• 

[17] Krnglov, S. 1. 2014, Phys. Rev. D, 89, 064004 (arXiv:1310.6915). 

[18] Liddle, A. R., and Lyth, D. H., 2000, ”Gosmological Inflation and Large- 
scale Strncture”, Gambridge University Press: Gambridge 

[19] Magnano, G., and Sokolowski, L. M. 1994, Phys. Rev. D, 50, 5039 
(arXiv:gr-qc/9312008). 


18 


[20] Muller, V., Schmidt, H.-J., and Starobinsky, A. A. 1998. Phys. Lett. B, 
202, 198 

[21] Nojiri, S., and Odintsov, S. D. 2003, Phys. Rev. D 68, 123512 
(arXiv:hep-th/0307288). 

[22] Nojiri, S., and Odintsov, S. D. 2007a, Phys. Lett. B 657, 238 
(arXiv:0707.1941). 

[23] Nojiri, S., and Odintsov, S. D. 2007b, Int. J. Geom. Meth. Mod. Phys. 
4 115 (arXiv:hep-th/0601213). 

[24] Nojiri, S., and Odintsov, S. D. 2008, Phys. Rev. D, 77, 026007 
(arXiv:0710.1738). 

[25] Nojiri, S., and Odintsov, S. D. 2011, Phys. Rep., 505, 59 
(arXiv:1011.0544). 

[26] Nojiri, S., and Odintsov, S. D. 2015, Astrophys. Space Sci. 357, 39 
(arXiv:1412.2518). 

[27] Starobinsky, A. A. 2007, JETP Lett., 86, 157 (arXiv:0706.2041). 

[28] Starobinsky, A. A. 1980, Phys. Lett. B, 91, 99 


19 


